Abstract 



To reveal how nonequilibrium physics and relativity theory intertwine, 
this articles studies relativistic Brownian motion under cosmic expansion. 
Two fluctuation theorems for the entropy As, which is locally produced 
in this extreme nonequilibrium situation, are presented and proven. The 
flrst, ^e~^°) = 1, is a generalization of the second law of thermody- 
namics, that remains valid at relativistic particle energies and under high 
cosmic expansion rates. From this relation follows, that the probability 
to observe a local reduction of entropy is exponentially small even if the 
universe was to recoUapse. For the special case of the Einstein-de Sitter 
universe an additional relation, i^e"^""^'') = 1, is derived which holds 
simultaneously with the first relation and where A/i is proportional to the 
Hubble constant. Furthermore, the fluctuation theorems are shown to 
provide a physical criterion to resolve the known discretization dilemma 
arising in special-relativistic Brownian motion. Explicit examples and a 
general method for the computation of non-Gaussian entropy fluctuations 
are provided. 
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1 Introduction 

The physical basis of the direction of time has been discussed since Boltzmann's 
H- Theorem in 1872. A priori, the thermodynamic arrow of time has to be 
distinguished from the possibihty of a prime direction of time defined by the 
expansion of the universe jl] . By now we know that due to the dominating dark 
energy component of about 72%, our universe is very hkely to expand forever 
[2]. Yet, the fascinating cosmological arrow of time could not be based on firm 
theoretical ground [H IH [5l [6l [7] . So one may still ask: Is it a mere coincidence 
that our memory strictly refers to times when the size of the universe was 
smaller? Put in physical terms, the guiding question of this article is: Does 
the cosmic expansion rate effect the production of entropy by nonequilibrium 
processes? While a general theory of nonequilibrium thermodynamics does not 
exist, fluctuation theorems (FTs) provide a unique starting point to develop the 
means to address such a fundamental question. First progress in the description 
of entropy production and giant fluctuations beyond linear-response was made in 
the 1970s j8j , and major advance was achieved in recent years with the derivation 
of FTs for various classes of systems UHl HH [HI [131 HI] • FTs generalize the 
second law of thermodynamics. The second law states that the Gibbs entropy 
S of an ensemble may not decrease, 

AS" > at any time. (1.1) 

The FTs naturally extend the concept of entropy and allow statements about 
the probability to observe isolated "violations" of (|1.1|) . To this end the FTs 
assign a change As of entropy to an observation of few or even single particles. 
When a nonequilibrium system of finite size is observed, the entropy As pro- 
duced within a certain time interval is a fluctuating quantity. The founders of 
statistical mechanics, L. Boltzmann and J.W. Gibbs were well aware that the 
second law holds only for the entropy AS" = (As) of an infinite ensemble. The 
angle brackets denote the ensemble average over observations of equal systems. 
Boltzmann mentioned the possibility of "violations" , As < 0, in his famous 
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reply to the Poincare recurrence objection (in a written argument with E. Zer- 
melo) and designated the second law as a theorem of probability ( "Wahrschein- 
lichkeitssatz" ) , emphasizing that the second law cannot be expected to hold for 
few particles US]. Loco citato, he referred to Gibbs [TB] who had concluded: 
"The impossibility of an incompensated decrease of entropy seems to be reduced 
to an improbability." It is this improbability that is quantified by FTs. 
For the steady state of strongly chaotic systems the detailed FT, 



was proven in the limit of infinite observation time [10] . The detailed FT (|1.2p 
was also derived in [11] for a non-relativistic particle in contact with a heat 
bath at temperature T. Initially, only the external change in the bath entropy, 
Ase = AQ/T with the energy AQ dissipated into the surrounding bath, was 
taken into account [TTJ [T^] . In [T3] it was pointed out that when the particle 
is assigned an intrinsic entropy Sg = — /cb InP (with the particle's phase space 
density P), the sum of intrinsic and external entropy. As — Asg + Asg, obeys 
the FT (|1.2p even for finite observation tim^ll. This is the definition of entropy 
applied throughout the article. Furthermore, for non-stationary states in the 
presence of time-dependent driving forces, an integral FT of the form 



was proven and linked to the Jarzynski relation (141 . Technically, the averaging 
over observations, {■■■), is a path integral over trajectories, which is explained 
in the following section. 

The detailed FT (fL2| . the integral FT pTS]) . and the second law of thermo- 
dynamics (|l.ip form a consistent hierarchy of statements: from Eq. (|1.2p follows 
(|1.3p by integrating over a, and Eq. (|1.3p implies (jl.ip by virtue of the Jensen 
inequality. 

We proceed as follows. In Sec. |2j FTs for general processes are derived. 
The recent unification [TT] [181 [IS] of Einstein's 1905 publications on Brownian 
motion [20] and special relativity [211 is briefly reviewed in Sec. [3] Based on 
these findings, the results of Seed] are twofold. First, for the relativistic Brow- 
nian processes of [TH] and [Hj , we reconcile the FTs (|1.2p and (|1.3p , which have 
become a paradigm of nonequilibrium physics, with special relativity. For the 
similar process suggested in [17], FTs follow by analogous reasoning. In [18] 
and [19] it was pointed out that the relativistic time dilation leads to multi- 
plicative coupling, necessitating a careful choice of the discretization rule. We 
show explicitly that there is one relativistic detailed FT (|1.2p and one relativis- 
tic integral FT (|1.3p valid for all choices. Second, we shall find the physically 

^ In the earlier FT of Evans and Searles a similar term was added to the dissipation function 
(Eq. (2.6) in [22]), that is not present in the Gallavotti-Cohen FT [10], resulting in a FT for 
chaotic systems which holds for finite time. For the stochastic formulation of finite time FTs 
such a term was considered in [23], Eqs. (5.9) and (5.11). 



Prob(As = +a ks) 



e" for any a. 



(1.2) 



Prob(A.s = —a /cb) 




(1.3) 
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correct expression for the entropy production following from relativistic FTs 
when the Hanggi-Klimontovich discretization rule is applied. In Sec. [5] we go 
beyond special relativity with a set of two general-relativistic integral FTs for 
the cosmological standard model. These expose clearly the role of cosmic ex- 
pansion in entropy production. We shall identify the entropy production which 
is solely due to the Hubble expansion of space. Such entropy producing pro- 
cesses dominate when the expansion rate of the universe exceeds the particle 
scattering rate, for instance in an early inflationary phase after the big bang. 
The Secs.|3]and[n]conclude each with examples where we explicitly compute the 
non-Gaussian fluctuations P(As), P(Ass) and P(Asc) of entropy production. 

2 Stochastic formulation of the fluctuation the- 
orem 

This section gives a general derivation of the FT for stochastic processes and 
emphasizes that every broken symmetry implies a FT. 

2.1 The integral fluctuation theorem 

Let T{t) denote the state of the observed system, which performs a time con- 
tinuous stochastic process under the influence of a thermal environment. Its 
stochastic dynamics are described completely by the probability distribution 
P[r,C], which gives the probability to observe a certain system trajectory T. 
This probability depends on the environmental conditions: C describes a set 
of external parameters, such as the environmental temperature T{t), external 
forces (for example acting on a charged system by an electric held E(<)), or 
- as we shall consider finally in the general-relativistic case ~ the curvature of 
spacetime. All these external parameters may vary during the process, so that 
C{t) is a deterministic protocol. The general idea underlying stochastic formula- 
tions of fluctuation theorems is as follows. Consider an arbitrary transformation 
T, which docs not leave the physical dynamics (represented by a Langevin or 
Fokker-Planck equation) invariant. While P[r, C] describes the dynamics of the 
original stochastic system, the transformed stochastic dynamics will be given by 
another probability distribution P[r, C]. Assuming P and P to have the same 
support, we define 



to quantify the symmetry breaking of the transformation T for every trajectory 
r. For the quantity As defined in (|2.ip . an integral FT of the form (|1.3p is a 
mathematical identity, 



As[r,C] EE fcBlnP[r,C] - fcBlnP[r,C] , 



(2.1) 




The path integration J T) [F] covers all continuous functions F, weighted by the 
probability P[F,C]. 
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We are interested in a fluctuation theorem that quantifies the irreversibility 
of the stochastic process. This is why we choose the transformation T to be 
time reversajl, 

p[r,c] -p[f,c], 

with f (+<) = r(-i) 

and C(+<) C(-t) for all t. 

To get a result on the total entropy production, the transformation T acts 
globally by reversing both, the stochastic system trajectory T and the time 
dependence of the environment C. 

The probability P[r, C] to observe a stochastic trajectory [Fjl^J! in the time 
interval (— r, +t) depends on the initial conditions, which are given by r(— r) 
for a Markov process or by the history [r]Z^_/p for a system with memory time 
T (which may be infinite): 

P[r, C] = Pi„ Pp , (2.2) 



with Pi, 



P(r,C)|_^, ifMarkovian 
P[r, C]I^_y, if with memory 

and Pp = P([r,C]t; |in) . 



We refer to the initial state or the history as the in-state of the system, which 
is distributed according to the first factor Pin in (|2.2p . The second factor Pp is 
the (forward) propagator on the time interval (— t, +r) under the influence of 
the thermal environment. Analogously, the time- reversed probability is written 
as 

P[r, C] - P[f , C] - Pout Pk ■ (2.3) 
Inserting (|2.2p and (|2.3p in (|2.ip . As decomposes into the sum 

As = Ass + Ase , (2.4) 

with Ass = -fee hi Pout + In Pin, (2.5) 
Pp 

and Asc — fee In — . (2-6) 

The first term (|2.5|) is the change of the system entropy Ss = —fee In P, as the 
system state changes from "j = in" to "j = out". The expression Ss(Pj) = 
— fceluPj for the system entropy was suggested in [13] for a Markov process 
and is a widely accepted definition because Sg resembles the Boltzmann entropy 
and Ss = (sg) = —fee (InPj) coincides with the Gibbs entropy of the ensemble 
Pr 

The second term (I2.6P is a Crooks relation [53] defined by the forward and 
reversed time evolution under the stochastic influence of the thermal environ- 
ment. We have to show that Aso as introduced in (|2.6p equals exactly the 



^If r is the phase space vector (x, p), the momenta are inverted, T{t) = (x(— t), — p(— t)). 
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entropy produced in the thermal environment, so that As = Asg + Ase is the 
total entropy production. It is the objective of this article to evaluate (|2.6p 
for a thermal environment at relativistic energies to adjudicate on the physical 
interpretation as environmental entropy. 

For Markov processes, such as the relativistic Brownian motion discussed in 
the next section, the forward and reverse propagators are infinite products of 
transition probabilities. 



Pf = lim J]nta,';.(rfc_iK.r,.,C(Aifc)) and 

n 

Pr = lim []Pttt(rfc^rfc_i,C(Atfc)) , 



n — >oo 

k=l 



so that entropy production is local in time: For the environmental entropy 
follows 



Aso = / Sc{t) dt with 



pdlns(f-^f+,C(t)) 



Se(t) dt = fcB In ^T^TTTTT ' (2-7) 



and the change in system entropy is Asg = Ss(+t) — Ss(— r) with Ss{t) = 
-kB lnP(r(t),C(i)). The probability density P(r, i) = P{T,C{t)) evolves 
according to the continuity equation, 

atP(r,i) + Vroj(r,i) = , (2.8) 

with the probability current 

j(r,i) = ^^ — p^oM„+i(r,i) P(r,t) . (2.9) 

n=a 

The Helfand moments M„ are tensors of order n which are related to the tran- 
sition probability P^^i^^^ by [55] 

M„(r, t)dt^J (r' - r)" p,f_ (r ^ r', c(t)) dr' . 

The higher moments are present only if the heat bath in which the system is 
embedded is out of equilibrium. A possible system for relativistic Brownian 
motion is an electron which couples by Compton scattering to a gas of photons. 
We assume that such a heat bath is in local equilibrium so that we have a 
well-defined temperature T(x, <) yielding an isotropic diffusion M2 oc TI with 
vanishing higher moments, M„ = for n > 2. Equation (|2.8p then reduces to 
the Fokker-Planck equation and the transition probabilities are Gaussian. 
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2.2 The detailed fluctuation theorem 



While the integral FT derived above holds for arbitrary environmental condi- 
tions C, the stronger detailed FT (|1.2p holds if the deterministic protocol is 
invariant under time-reversal, C = C. The general derivation for the quantity 
As defined in (|2.ip is also done conveniently by path integration. The proba- 
bility to observe a production of entropy As — a is 

Prob(As = a fee) 
= J P[r, C] S{As[r, C]=a ks) V[T] 

= J P[f, C] e'^'/''^ S{As[T, C] = a fee) V[T] 

= e'^ y" P[f, C] 6{As[T, C]=a fee) 2?[r] 

= e'^ y P[f, C] (5(As[f , C] = -a fee) 2?[r] . 

In the second and last equality we exploited Eq. (|2.ip . Using the trivial fact 
that the path integration can be reordered in time, we arrive at 

Prob(As = a A;b) = e" J P[T, C] 6{As[T, C] = -a ks) V[T] . 

Comparing this result with the probability 

Prob(As = -a kn) = J P[T, C] S{As[T, C] = -a fee) V[T] 



to observe a reduction As = —a ks, yields the detailed FT (|1.2[) for any sym- 
metric protocol, C = C. Therefore the detailed FT holds not only in the steady 
state, which the system reaches under time-independent forcing, C{t) — const, 
but also for example in periodically changing conditions that are symmetric 
with respect to the observed time-frame {—t,+t). 

We conclude this general derivation of FTs with the remark that the pre- 
sented formulation gives a unifying perspective on the distinct FTs of [13 and 
pS] , Equation (6) in |2B] is generalized by Eq. (|2.ip . while the Eqs. (|2.5p and 
(|2.6p correspond to the decomposition of entropy according to the Eqs. (5) and 
(14) in [13j respectively. 



3 Relativistic Brownian motion 

The derivation of the FTs (|1.2p and (|1.3p in Sec. [2] uses the abstract expression 
(|2.7p for the entropy production in the embedding heat bath. As emphasized 
before, this expression has to be evaluated for a physical process to allow for a 
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physical interpretation as entropy. An instructive process is relativistic Brown- 
ian motion. 

To minimize technicalities, we consider first the one- dimensional special- 
relativistic motion of a particle with rest mass m in a heat bath at temperature 
T. The generalization to higher spatial dimensions is straightforward. Even 
if we would allow the particle to equilibrate with its environment, the mean 
squared velocity may not obey the non-relativistic law (v^) = k^T /m in the high 
temperature limit, since the finite speed of light defines an insurmountable upper 
bound. The special-relativistic nonequilibrium Brownian motion, giving rise to 
bounded velocity distributions, has been set forth in [T71[T51[Tn] using both, the 
language of stochastic differential equations (relativistic Langevin equations) 
and the language of probability densities (relativistic Fokker-Planck equations). 
Simulations of this relativistic stochastic process have been applied to analyze 
scattering experiments of quark-gluon plasma [27j . As in the familiar non- 
relativistic case |1S], a deterministic force acts on the particle in the rest 
frame of the heat bath, 

dpd = Fddt = -vp dt , (3.1) 

so that the time scale of dissipation is In the relativistic generalization (|3.ip , 
the non-relativistic momentum mv is replaced hy p = = mv j ^\ — jc? ^ 
which is the spatial component of the relativistic momentum vector p". As 
common, Greek indices refer to temporal (a — 0) and spatial components. The 
signature of the Minkowski metric tensor is r\a^ = "rf^ = diag(— 1, 1). Moreover, 
Einstein's summation convention is invoked throughout. Since the rest mass is 
not altered in elastic collisions, p°'pa — —{mc)'^ ~ const, the change in the 
momentum vector dp" is always "orthogonal" to Pa in the sense of 

Pa^dp" = . (3.2) 

This means that the classical particle cannot leave its mass shell p'^Pa = — (jnc)'^ , 
which is nothing but its dispersion relation, 

E^p'^c^^imc^f + ipc)^ . (3.3) 

The general solution of p.2p is the projection dp" — {S^ +p"pi3/{mc)'^)^^ of an 
arbitrary Lorentz vector It is readily confirmed that the choice 

dPd = -mz^ + ^) Chdr (3.4) 

reduces to Eq. (13. ip in the rest frame of the bath with the bath velocity vector 
''^bath — i^' ^) particle's proper time r. Hence, Eq. (|3.4P is the generalized 

Lorentz-invariant deterministic part of the Brownian motior|j. 

The description of relativistic Brownian motion is completed by Lorentz- 
invariant stochastic changes dp" of the momentum caused by the impacts of 

•^Equation (16) in 18 contains an identically vanishing term. 
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the surrounding heat bath at temperature T. The derivation is guided by two 
principles: first, the relativistic momentum is the proper quantity performing a 
Wiener process, since it is physically exchanged and additive, whereas the veloc- 
ity is well-known not to be additive in special relativity. The second postulate 
demands that the distribution is Gaussian in the instantaneous rest frame of the 
particle. This connects the relativistic Brownian motion to the non-relativistic 
case. These principles determine the exchanged momenta dp" to be distributed 
according to (cf. Eq. (35c) in [TH]) 

Pcou b'^,dpj = — ^=^exp I — — 1. 3.5) 

The Dirac distribution 5{pj3dp^) in (|3.5p guarantees that the mass-shell con- 
dition (|3.2p is also fulfilled by the stochastic impacts, since they are elastic. 
While the relativistic momentum p is additive and unbounded, the velocity is 
restricted to the open interval (— c, -\-c). This can be seen by the elegant relation 
V /c^ = p/ E in the rest frame of the bath, which is equivalent to 

Ax = dt . (3.6) 

^ [mcY + p"^ 

As mentioned before in the context of the general Kramers-Moyal expansion 
(|2.9p . the bath temperature T is defined by the Einstein relation, 



^ = kBTmv , (3.7) 
with the momentum diffusion constant !^ (cf. Eq. (59) in ,18j). 



4 Relativistic fluctuation theorem 

We have now the manifestly Lorcntz-invariant Langevin equation 

dp'^^dp^+dp^ (4.1) 



with the deterministic part given by p.4p and the stochastic part described by 
(|3.5p at hand. Specifying (|4.ip to the rest frame of the bath yields 

dp = -up dt + dps . (4.2) 

The probability density of the exchanged momenta dps is found by integrating 
out the dpg -component in (|3.5p . cf. [T5| : 



exp(-dp3V(4i^y'l+(;£^dt) 
Pcoii(p, dps) = ^ , ■ (4-3) 

This exhibits the discretization dilemma: A discretization rule has to be imposed 
on (j4.3p since relativistic invariance does not determine whether p in (j4.3p refers 
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to the particle momentum p_ before the colhsion (pre-point rule of Ito), to the 
post-point p+ = p- -\-dp (Hanggi-Klimontovich), or to the midpoint (p_ +p^)/2 
(Fisk-Stratonovich) . 

The Eqs. p.6p , (|4.2p and (|4.3p establish the relativistic stochastic motion of 
the Brownian particle in phase space. The corresponding transition probability 
is uniquely determined by the discretization rule: 

/.^. + d. \_ <d.x-^d^ 



P' 

trans 



p>-^p + dp J 2y^TrSiEdt/mc^ 



dp + updt^{l^K)-^^dt 



2 



I \ rnc^ dp , , 

" A^Edt/mc^ ) ■ (^-^^ 

The discretization is contained in the parameter k, < k < 1. Hanggi- 
Klimontovich, Fisk-Stratonovich, or Ito correspond to the values k = 0, ^, 
or 1 respectively. 

Let us now investigate the consequences for entropy production arising out 
of the special-relativistic discretization dilemma. As derived in Sec. [21 the total 
entropy is a sum of the particle intrinsic entropy Ss = —fee In P with the parti- 
cle's nonequilibrium phase space density P{x,p, t), and the external entropy So 
of the ambient heat bath at temperature T. 

Inserting the probability current (|2.9p in momentum space 



■ / .^ f ^ , , &EdP{x,p,t) 

jp{x,p,t) = - lyp + K -— P{x,p,t) yil^ (4.5) 

\ TOC^ dp / mc^ op 

in the differential dsg of the particle entropy Ss{t) — — /cb In P(a;(t),p(t), t) we 
find the equation of motion (generalizing Eq. (7) in [13J) for Sg, 

dss = dssU=o + K fcs dlni; . (4.6) 

Here we have isolated the second term which depends on the discretization rule 
applied. 

The entropy production dsc in the bath follows by contrasting the transition 
probabilities of the trajectory F — {x,p) with its time-reverse T — {x,—p) to 
extract the irreversible part, InP^l^^ (F" i-^ r+) -lnPtf^„g(f " i-^ f+), causing 
the dissipation p.7p . From a brief computation we find: 



dsp — fcn In 



trans 



X 1-^ X + dx 
p 1-^ p + dp 



pdt 
trans 



X + dx X 
—p — dp i—i^ —p 



d E 

= -—-nkBdluE. (4.7) 

The Eqs. (|4.6p and (|4.7p reveal that although the relativistic Brownian mo- 
tion is physically incquivalent depending on k, the fluctuations of the total 
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entropy s = Sg + So are independent of k. Explicitly, the change of the total 
entropy is 

ds d\nP d\nP mc^jp 

— = t: — dt da; + ^^„ dv . 

kB dt dx 2>EP ^ 

Two technical comments are here in order. First, when computing dse in 
(|4.7p the notation has to carefully distinguish between initial r~ and finale state 
r+, and one should consider the quotient F(r" r-+dr)/P(r~+dr T") 
as done in (j4.7p . Writing the back transition in the numerator in the form 
P(r+ r+ — dr) would be correct yet unfavorable for evaluation, because 
common F^-factors could not be cancelled out. In transforming P(r+ r+ — 
dr) to P(r^+dr r~) the known spurious drift of the multiplicative coupling 
has to be taken into account [53]. Second, the discretization term in (|4.7p can 
be absorbed by defining a more complicated fluctuation-dissipation theorem, 
however in this article we use exclusively the Einstein relation (|3.7p . 

The path integration of the results (|4.6p and (|4.7p according to Sec. [5] yields 
the detailed FT (|1.2p for time-symmetric environments, and the integral FT 
(jl.3p for arbitrary environmental conditions, with entropy fluctuations As ob- 
served over finite time. Therewith we have proven relativistic FTs that are 
unaffected by the discretization dilemma. 

Furthermore, we are now in a position to address the physical choice of k, by 
virtue of the FT. Because of energy conservation, the energy —dE in (14. 7|) lost 
by the particle equals the heat dQ gained by the ambient bath: 

dse ^ - K fee dln£; . (4.8) 

In the non-relativistic regime the particle energy E = mc^ + -Ekin is dominated 
by the energy of the rest mass m so that the second term in (|4.8p vanishes for 

d In £; = — 1^ ^kin , 

and we recover the non-relativistic FTs [13] ■ At arbitrary relativistic energies 
p.3p the Hanggi-Klimontovich rule, k = 0, entails the correct expression for the 
entropy 

d^e = ^ , (4.9) 

which is produced in the heat bath. 



4.1 Generalizations in the framework of special relativity 

To generalize the FTs to n spatial dimensions, momentum and force in Eq. (|3.ip 
are simply substituted by their spatial vectors and the Greek indices in the 
Lorentz-invariant Eqs. (|3.4p and (|3.5p take values up to n. After integrating 
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Figure 1: A Brownian particle commuting between different thermal environ- 
ments. The non-Gaussian fluctuations of entropy occurring in this system can 
be evaluated analytically (cf. Figs. [2] and [S]). 



out the temporal component p", the distribution (|4.3p is found to contain a 
quadratic form A instead of the square in the exponent (cf. Eq. (15) in [TOj ) 
with tensor components 

A,j = % - -^P^Pj ■ (4.10) 

The FTs follow using the fact that p is an eigenvector of A. No complications 
are caused by allowing an inhomogeneous heat bath, where the temperature 
T and the dissipation rate ly vary in space. As far as the integral FT (|1.3p is 
concerned, a bath temperature evolving in time is also permitted (as part of the 
environmental condition C(t) in Sec. [5]). Since the time-asymmetric part enters 
(j4.9p . the dissipation rate v may be an even function of the momentum, z^(p) = 
z/(— p). This is of physical relevance since i/ is known not to be constant even 
for most non-relativistic processes |2H]- As mentioned in the general derivation 
of Sec. m an arbitrary time-dependent external force Fe{t) (being also part of 
the environmental condition C{t) defined in the rest frame of the bath) does 
not pose a problem. After adding Fe{t) to the deterministic force Fd in (|3.ip 
we find the expression dse = dQ/T with the heat dQ = —dE + Fedx. This is 
the first law of thermodynamics stated in the frame of the bath. 

4.2 The commuting Brownian particle 

We give (to the author's knowledge) the first example where the non-Gaussian 
fluctuations of particle entropy Asg, environmental entropy Aso, and total en- 
tropy As can be evaluated exactly. A complementary method which allows the 
general numerical computation of fluctuations by iteration will be proposed in 
Sec.lOl 

Consider two heat baths at temperatures Ti and T2 with a Brownian particle 
moving initially in Ti. After the equilibration time, its momentum pi (in units of 



12 



mc) will be distributed according to the Jiittner-MaxwcU distribution ip{pi,Ti) 
[291. where 



= (4.11) 

Z{T) ^ ' 

is the equilibrium solution of the Brownian motion presented in Sec. [31 The 
last formulation in (|4.1ip (following from a/1 + _ \ = p"^ yj\ is 
convenient for the low momentum limit. To keep formulas concise, we measure 
heat in units of mc^ and entropy in units of fee, so that temperature T is 
measured in units of mc^ /k^. The relativistic partition sum Z{T) equals 

z(r) = 2 exp(i/r)Xi(i/r) , (4.12) 

with Ki being the first modified Bcssel function of the second kind. The Brow- 
nian particle can pass to the bath T2 through an opening (cf. Fig. [T]) . This 
opening is small enough to keep the baths at different temperatures and to en- 
sure that the Brownian particle spends enough time in T2 before returning to 
Ti . So its momentum pf on return has become uncorrelated to the initial value 
Pi and is distributed according to tf{pf, T2). The change of the particle entropy 
Asg and the environmental entropy Aso during the relaxation of the Brownian 
particle in T2 can be expressed using tp: 

As, = s^iU) - s^iU) =\n^^p^ (4.13a) 



^{PhT2) 



As. =^ =1-^)7^ (4.13b) 



As =As, + Ase =ln4^il^ (4.13c) 

The total entropy As in (|4.13cp follows from the above definitions of particle 
entropy (|4.13ap and environmental entropjQ (j4.13b[l . From the resulting expres- 
sion (|4.13cp we find the macroscopic Gibbs entropy, 

AS = (As) = / ^(pi, Ti) In ^^f^ dpi , 

to equal the relative entropy of the baths, 

A5 = 5kl(Ti||T2) , (4.14) 
which is also known as the KuUback-Leibler distance ^301. 



^Aside from the physical expression AQ/T used in I l4.13bll . the expression 
In (<^(pf, r)/(/j(pi, T)) in terms of the equilibrium distribution ip is directly related to 

the definition 1)2. 7| l by the principle of detailed balance, ip(pi,T) P^^^^{pi pi,T) = 

'p{pi, T) -Ptrans(Pf Pi' -^)' because the baths themselves are in local equilibrium. 
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The trajectory entropies (|4.13p depend only on the pair {pi,Pf) of initial and 
end point in momentum space. Therefore the distributions P(As), P(Ass) and 
P(Asc) follow not from path integrals but ordinary integrals such as P(Asc) = 
/ (p{pi, Ti)(p{p{,T2) S (Asc - Asc(pi,Pf)) dpidpf, where the expression Ase(pi,Pf) 
(|4.13b[) is inserted in the Dirac delta function to sum over all trajectories yield- 
ing a certain entropy increment Ase- Because of its physical relevance, we begin 
with the explicit non-relativistic results, T <C mc^ /k-Q, when 93 (j4.1ip becomes 
the Maxwell-Boltzmann distribution: 

P(A.) = ^M^l^,-^ (4.15a) 
^ttA (A,s - so) 

TT 

P(Ase) = ^e^^-^^ ifo(^|Ase|^^) (4.15c) 

The abbreviations A ~ a^^ — 1 and so = ^Ina contain the dependence on 
the temperature ratio a — T2/T1. The Heaviside step function is denoted by O. 
The distribution functions (|4.15p are plotted in Fig. [21 With the Bessel function 
Kq appearing in (|4.15p . the distributions for Asg and Ase have logarithmic 
divergences at sq and respectively. The distribution of As has the stronger 
inverse square root divergence as As approaches sq from above and vanishes 
below Sq. From (|4.15ap the integral FT (|1.3p can be verified directly, while the 
detailed FT (|1.2p is obviously not fulfilled (as it has to be since the embedding 
temperature for the Brownian particle changes randomly with time) . We remark 
that (I4.15ap is not simply the convolution of (I4.15bp and (|4.15cp because Asg 
and Asc arc highly correlated. 

The macroscopic entropies ASs = (Asg), AS'c = (Asc), and AS* = (As) = 
ASs + ASc > are the mean values of the distributions (|4.15p : 

In a , , 

AS, = so = — (4.16a) 

A^e = (4.16b) 

After the Brownian particle has visited both reservoirs once, the total macro- 
scopic entropy increment has the symmetric form 

AS{Ti ^ Ta ^ Ti) = AS{T2 ^ Ti ^ T2) = > . (4.17) 

In the relativistic regime, mc^ defines a third energy scale, so that the results 
no longer depend only on the ratio of temperatures. The singularity at As = sq 
is shifted to the position 

so^mfj^ (4.18) 
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Figure 2: The exact expressions (|4.15[) for the distribution of particle entropy 
Ass (blue), environmental entropy Aso (green), and total entropy As — Ass + 
Ase (red). The plot is for Ti = AT2 ^ mc^/kB- The inset shows the dependence 
on the temperature ratio for the macroscopic entropies ()4.16|) : ASs — (Ass), 
AS'e = (Aso), and AS* = A5s + ASc which is non- negative according to the 
second law of thermodynamics (jl.ip . 



in terms of the partition sum (j4.12ll . The first relativistic correction of the 
partition sum is 



Z{T) = V2^ ^1 + + 0{T'^)^ 



(4.19) 



so that So depends on the temperature difference AT = T2 — Ti in first order: 
so-iln^ + ^Ar + 0(T2,r|) . (4.20) 

In the ultra-relativistic regime, T 3> mc^/fce, the partition sum becomes linear 
in T, 

Z{T) = 2T + 2 + 0{l/T) , (4.21) 

so that the position sq of the singularity depends on the ratio of temperatures 
as in the non-relativistic limit and reaches twice its non-rclativistic value, 

T 

lim So = 2 lim so = In — . (4.22) 

The relativistic distribution functions are sums of Bessel functions. For exam- 
ple the system entropy As is distributed at arbitrary temperatures Ti and T2 
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according to 

The normalization factor in (|4.23p is 

W,T.) = ^g^ (4.24) 

and the function / in (|4.23p is defined by the integral 

. ^ -z r ^ l + ax/2 l + b{x/2 + z) 

f(a,b,z) = e ^ / da; ^ , ^ — , ^ . (4.25) 

Jo Vi\/^T2l ^l + aa;/4 y^l + b{x/4: + z/2) 

The non- relativistic limit (|4.15bp follows from /(O, 0,2) = Koiz) and /^(0,0) = 
TT. The first relativistic corrections are 

P{As) = _ goiAss-so) 
+ T 5i(Ass-so) 

+ AT 52(As, -so) +0(T2,T|,riT2) , (4.26) 

with the mean temperature T = (Ti + T2)/2 and the temperature difference 
AT = T2-T2. We remark that the functions 

9oiz) = ^ (4.27a) 

TT 

.gi(z) = |-(|z|ifi(|z|)-ifo(|2|)) (4.27b) 

52(^) = ^^ifo(kl) (4.27c) 

exhibit the symmetry Ti T2 of the system, gi{—z) = gi{z), g2{—z) = —52(2:), 
and preserve the normalization at any order, gj{z)dz ~ Sqj- 

In the ultra-relativistic limit, we find N 2y/T{I^ ^^M^ and /(Ti, T2, z) 



y/TYT^c-'' g^Sl), so that 



g-|ASa-So 

P(Ass) = (4.28) 



is an exponential distribution. It is only in the ultra-relativistic limit that 
the logarithmic divergence at sq vanishes in favor of a kink (cf . Fig [3]) . In 
the intermediate relativistic regime (fceT « m(?) the distribution P(Ass) has 
skewness. The exact distribution (I4.23P is shown for a fixed temperature ratio 
Ti = 4r2 as the geometric mean \JT\T2 is increased form zero (non-relativistic 
limit) to infinity (ultra-relativistic limit) in Fig [31 
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Entropy / fee 

Figure 3: This plot shows the distribution of system entropy, P(Ass), as we pass 
from the non-relativistic regime to the ultra-relativistic regime. The parameter 
attached to each graph is fee -sJTiTij (mt?) which assumes the values zero (non- 
relativistic limit), 1, 10 and oo (ultra-relativistic limit). As we approach the 
ultra-relativistic limit, the mean doubles and the spread of fluctuations widens, 
but does not diverge. The logarithmic peak reduces to a kink. Note that we are 
discussing the classical relativistic regime. Quantum corrections, depending on 
the particle spin, are expected when pair creation sets in. 
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Figure 4: A sketch of spacetime showing a spatial shce of the heat bath at fixed 
time and the world line of a Brownian particle in a (locally) expanding universe. 



5 Generalizations in the framework of general 
relativity 

The monotonic increase of entropy is a fundamental principle of physics and 
the universe is known to expand, as was discovered by E. Hubble in 1929. 
The discussion whether there is a direct connection between these observa- 
tions has never stopped [31 21 O [7]. Therefore we aspire a formulation 
of the FT consistent with general relativity, but we restrict ourselves to the 
class of Friedmann-Lemaitre models, which describe a spatially homogenous 
and isotropic, expanding or contracting universe. The corresponding line ele- 
ment (given by the Robertson- Walker metric) is —dt^ + dr^. The important 
difference compared to special relativity is that the spatial part, dr^, is scaled 
by a time dependent factor R(t) describing the expansion or contraction of the 
universe: 

dr^ = R^t) dC d^^ . (5.1) 

The Latin indices describe spatial components numbered by 1 to 3. We do 
not have to deal with the details of the metric tensor h describing the spatial 
curvature. The result will be valid for all possible geometries. The expansion 
rate H{t) — R{t)/R{t), named Hubble function, is one of the most important 
quantities in cosmology and its present value is a direct observable [3T]. The 
typical frame for a cosmic heat bath is the frame of the cosmic microwave 
background. 

5.1 Cosmological fluctuation theorem 

In general relativity, the correct equations of motion include the covariant dif- 
ferential Dp of the momentum. (Denoting by p the 4-vector, the components of 
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Dp are Dp" — dp" + r^^p'^da:".) Its spatial components replace the left hand 
side of (|4.2p and can be split up into a spatially covariant part, ('^^Dp, and a 
contribution due to the time-dependent scaling: 

Bp ^ '•^^Bp + H{t) p dt . (5.2) 

Therefore the covariant Langevin equation, generalizing Eq. (|4.2p to be valid in 
an expanding or contracting universe of arbitrary spatial geometry, reads 

(3)Dp = - [i^{p, t) + H{t)] pdt + (3)Dp^ _ (5.3) 

Herein, H enters as an additional damping term, which has caused the cooling 
during the expansion of our universe and is responsible for the cosmological red 
shift. The distribution of the stochastic impacts ^'^^Dps is found after substitut- 
ing hij for the Euclidean metric 5ij in (j4.10p . Applying the time-reversal map, 
we find that Eq. (|4.9[) gains a second term due to the cosmic expansion: 

dsc = — 7^ — It dlni? 
T ET 

The numerator (p, dr) in (|5.4p is the canonical line integral (canonical one- 
form) in phase space. The integral FT (|1.3p extends to an expanding {H > 0) 
or contracting (H < 0) spacetime when this second term is taken into account. 
It has a clear geometric interpretation: the Hubble function is the external 
curvature of space, 

BN = H dr , (5.5) 

with TV being the time-like normal vector to the space of the heat bath as 
depicted in Fig. [H This permits the second term in (|5.4p to be written as 

J ^{cosmic) 

{p,BN) 

Since the particle energy E = — — po — — (p, N) is the zero component of the 
4- vector p, the first term in (|5.4p equals the differential 



J (particle) 

_ d(p,iV) _ (Dp,7V) + (p,D7V) 

'^*e ~ T ^ T ' 

such that the sum of both terms is 

... ^ ^ (5^6) 

It is natural to think of the numerator (Dp, N) as the heat dQ ~ Tdsg exchanged 
with the bath, since it is the projection of the exchanged 4-momentum Dp on 
the local energy component N of the heat bath. 
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Cosmology is an example for the breaking of the first law, —dE = d{p, N) ^ 
(Dp, TV) = dQ, by non-static metrics. So we find ourselves in a remarkable situ- 
ation: There is no first law in cosmology, while the second law and furthermore 
the integral FT hold. 

The isolated cosmological entropy term dso'^°^™"^^ in (|5.4p would indeed un- 
dergo a change of sign if the expansion turned into contraction. But in the 
entire bath entropy (|5.6p the geodesic flow N enters as a projection, which does 
not imply a change of sign if N was to contract. Eventually, a decreasing total 
entropy s ~ Ss + Se is always exponentially unlikely as expressed by the integral 
FT (fTgl) . 

5.2 Additional theorems for the Einstein-de Sitter uni- 
verse 

So far, we applied the time-reversal transformation to arrive at (e"^'') — 1. 
As emphasized at the outset of the general derivation in Sec. [21 we are free to 
choose any other transformation from the mathematical point of view. Then the 
function in the exponent will no longer equal the entropy As. For instance, if the 
system is invariant under the chosen transformation, we will get the trivial result 
(e°) = 1. But for physically sensible transformations, the FT will remunerate 
us with non-trivial relations. In order to derive a FT that contains the cosmic 
expansion rate H, let us choose a local time-reversal transformation, which acts 
only on the local particle dynamics and leaves the sign of the global cosmic 
expansion rate H unchanged. Repeating the computation of Sec. 15.11 with the 
transformation H = H (local time reversal) instead oi H = —H (global time 
reversal) yields 

^g-(A.+Ah)^ = 1, with Ah = (5.7) 

The additional term Ah is proportional to the Hubble constant, the inverse 
temperature and the action A = J ((p, v) — AE) dt of the energy change AE = 
E/v. 

This demonstrates that the FT is an efficient technique to design relations 
that include those physical observables, which are most interesting for a given 
system or experiment. The second general relativistic FT (|5.7|) holds in addition 
to ()1.3p . At first glance one might be surprised that there is an infinity of FTs, all 
constraining the fluctuations of As. But since the distribution function P(As) 
is a point in the infinite dimensional (Banach) space of integrable functions, 
there has to be an infinity of physical constraints to determine P(As) uniquely. 

In many interesting stages of the cosmic evolution, such as the early (hypo- 
thetical) inflationary phase and the future phase of accelerated expansion, the 
size of the universe grows exponentially with time so that H is constant. During 
these periods the cosmic impact on the local relativistic Brownian motion with 
(|5.2p is time- independent. We can therefore immediately infer from the general 
derivation in Sec. l2.2[ that for these phases of the cosmic evolution the stronger 
detailed formulations of the FT hold as well. 
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5.3 The expanding universe 

The cosmological FTs of the Sees. 15.11 and 15.21 restrict the entropy fluctuations 
As caused by a relativistic particle. In this section we compute the detailed 
distribution of fluctuations explicitly for evolving cosmic environments. 

The entropy change ds ~ dsg + dsi'''*'^*"^''^' + dso'^°''™"^^ has contributions of 
the system entropy, Sg = — InP, and by heat exchange, dso = dQ/T. Since 
there is no first law, dE + dQ ^ 0, for the time dependent cosmic metric, the 
heat contribution dso splits up in a term due to the change of particle energy, 
^^(particie) _ jrj^ ^ .^^^^ ^ cosmological term, dso'^°™"^^ = —p^Hdt/{ET), as 
derived in Eq. (|5.4p . The method to compute distributions of fluctuations will 
be presented for the particle term, ^(p''''*"^''^) ^ which we abbreviate by Sp. It is 
straight forward to apply the method to the other terms. 

To compute the distribution P(Asp, Ai) of produced entropy Asp, we have 
to sum up dsp = —dE/T(t) over the observation time At. Therefore we have 
to evolve the process p{t) while book keeping the change of entropy Sp. This is 
done by extending the Fokker-Planck equation to evolve the joint distribution 
P{p, Sp,t). The evolution of entropy Sp is directly related to the dynamical 
variable p by the differential dsp — —dE{p)/T(t), since relativistic Brownian 
motion is restricted to the mass-shell (13. 3p . The similar evolution of P{p,E,t) 
is easily determined. One method is to include the Helfand moments (dE), 
(di?^), and {dE dp) into the probability current (|2.9p . from which wc find the 
Fokker-Planck equation dtP + dpjp + OeJe — for P{p, E, t). The correlation 
{dE dp) is important because dE is not independent from dp on the mass-shell. 
Equivalently, we can proceed using the Fokker-Planck equation dtP + dpjp = 
for P{p, t) with the current (|4.5p and substitute every differentiation dp by 
dp + so that the probability current is tangential to the mass-shell. After 

identifying ds^ = —TOe (j4.9p we arrive at the Fokker-Planck equation 

UtP ^ Fo{dp-^ds^)P 

= 'F„idp) - F,{dp) ds^ + F2 ■ dl] P (5.8) 
for the distribution P{p, Sp,t). The momentum operator is 
Foidp) = A + (A + T/E)pdp + ETdl . 
The entropic extensions of the Fokker-Planck Eq. (|5.8p are 

Fi{dp) = 1 + AP2 + 2pdp and F2 = p^/{ET) . 

The function A(t) = 1 -f H{t)/v contains the cosmic driving by expansion. This 
function of time is deterministic since we can safely neglect the back reaction of 
our tiny system on the cosmic evolution. The entropy fiuctuations P(Asp, At) 
follow from (|5.8p when solved for the initial condition 
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and after integrating out the momentum p: 

P{Asp, At) ^ [ P{p, Asp, At) dp . (5.10) 

JR 

The Fokker-Planck Eq. (|5.8p is solved by orthogonal functions. We expand 
the distribution P{p,Sp,t) in a series of Hermite polynomials with respect to 
the entropy dependence, so that the two-dimensional Fokker-Planck (|5.8p for 
P(p, Sp,i) reduces to an one-dimensional system for the coefficients ak{p,t). 
The coefficients ak{p, t) are simple linear combinations of the moments Mi{p, t), 



Mi{p,t) = J P{p,sp,t) 4 dsp 



so that the singular initial condition (|5.9p are represented by A/o(p, 0) = Po{p) 
and Ml {p, 0) = for alH > in a regular way. From (|5.8p follows after inte- 
grating by parts a hierarchy of differential equations for the moments Mi{p, t): 

^dtMi = Foidp)Mi + I Fiidp)Mi-i + l{l - 1) F2 Mi^2 ■ (5.11) 

The case / = reduces to the Fokker-Planck equation for the momentum, 
Mo{p,t) = P{p,t). Since (|5.1ip is a parabolic differential equation, numerical 
solutions for the Mi{p,t) can be obtained by standard techniques. Integrating 
p, we have the moments m;(Ai) = J Mi{p, At) dp for the distribution of en- 
tropy (|5.10|) . After computing iteratively a sufficient number of moments to;, 
the probability distribution for the entropy (|5.10p can be reconstructed by the 
algorithm presented in Appendix [A] 

Let us illustrate (|5.10p for a universe undergoing a transient inflation as 
sketched in the inset of Fig. \5[ Such a transition of the scale factor ranging from 
i?i to Rf according to 

, -Ri e-" + i?f e" 
R{t) = 

is a common toy-model for particle creation in quantum field theory [32j . The 
peak of the Hubble function shall be i/max, so that r = tH^^x/I- The inflation 
factor is / = 2{\/Rf — ^/Ri)/{y/R{ + \/R\). Neglecting quantum effects, the 
thermal heat bath, which may consist of photons or other massless particles, 
cools proportional to the inverse scale factor ^33^ , 

rp r) 

T{t) 



We choose the mean temperature Tmoan in the relativistic regime, feBTmcan = 
lOmc^. The universe inflates by the factor R{/Ri — 2. The cosmic forcing of 
the system depends on the ratio of the relaxation rate ly and the expansion 
rate ifmax- For the nonequilibrium distribution of the particle entropy Asp 
shown as solid line in Fig. [5] the dimensionless control parameter H„^ax/i^ equals 
100. As reference, the symmetric distribution of the relativistic equilibrium 
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Figure 5: The distribution of particle entropy Asp is shown for a cosmic in- 
flationary phase. The system starts in equilibrium at time t = —1/v with 
temperature Tj, undergoes a period of inflation centered at i = (cf. the inset), 
and equilibrates again until the time 1/v with the lower bath temperature Tf . 
The maximum of the Hubble function is -ffmax = lOOi^ for the distribution shown 
in solid line. The case of a static universe with zero mean entropy is plotted 
in dashed line. This system is described completely by the ratio -ffmax/ v of the 
cosmic expansion rate and the thermal relaxation rate, the temperature to mass 
ratio /cBT'moan/(™c^) = 10, and the inflation factor Rf/Ri = 2. 



with ffmax = is plotted in dashed line. When -ffmax/i^ assumes the values 1, 
10 and 100, the width Cs of the distribution P(Asp) increases monotonically, 
being equal to 1.45, 1.48 and 1.64 respectively. In contrast, the mean AS'p is 
not monotonic and assumes the values 0.68, 0.69 and 0.15 respectively. Five 
moments have been computed to construct Fig. [5l 



6 Conclusions 

Relativistic FTs have been established that remain valid for high temperatures 
or low masses, mc^ <^ k-oT. The integral FT, (e"'^'*) = 1, was found to hold 
also in the framework of general relativity as far as the cosmic expansion is 
concerned. 

With the additional FT ^e"'^*^^'') = 1 and the numerical example of 
Sec. 15.31 we can answer the question raised in the introduction: yes, the cos- 
mic expansion has an influence on the total entropy fluctuations As, and the 
mean values of individual terms such as the particle contribution Asi^'^'^*''^'''^ 
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Figure 6: A sketch of the space spanned by the functions (jA.3[) . The first n 
coefhcients 03, . . . , a„ are predetermine by the known moments (vertical hne). 
The higher coefficients a„+i, . . . are determined iteratively. The shaded region 
represents the non-Hnear (convex) space of non- negative functions. The inset 
shows the rapid convergence of the algorithm. 



can undergo a change of sign for cosmic contraction. However the relation 
(e"^") = 1 implies the second law, AS* = (As) > 0, so that for a macroscopic 
system the sign of AS" is independent of the cosmological evolution. 

On the theoretical road ahead, one may expect integral FTs to hold for ar- 
bitrary time-dependent and inhomogeneous fields, such as gravitational waves, 
when the concise expression ()5.6|) is applied. For the process originally intro- 
duced in [T7], the weaker inequality (|f .ip has been proven recently [31] under 
general conditions. 

Experimentally, the change of the environmental entropy Asc — ~AE/T 
can be measured by detecting single particles after a sequence of elastic colli- 
sions, i.e. collisions without decay or excitation of internal degrees of freedom. 
Such collisions are observed for heavy quarks (for instance the charm quark) 
which traverse the expanding quark-gluon plasma created by heavy-ion col- 
lisions. Nonequilibrium thermodynamical descriptions are common for these 
relativistic media [27]. The relativistic FT is not only subject of high energy 
physics and cosmology. The special-relativistic FT can be tested with a high- 
precision spectroscopy experiment by shining a laser on an excited granulate of 
glass or reflecting steal beads, so that the granulate serves as a heat bath and 
the photons are the ultra-relativistic "Brownian" particles. The environmental 
entropy Ase = —AE/T then follows from the measurement of the frequency 
shift Av = AE/h. 

The author is grateful for the comments of M. Brinkmann, J. Dunkcl, P. Hanggi, 
S. Herminghaus, C. Jarzynski, U. Seifert, and V. Zaburdaev. 
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A The truncated moment problem 



To reconstructing a distribution function, we are interested in an efficient algo- 
rithm that generates uniquely out of n > 2 given moments a continuous and 
non-negative function / on the real line, so that 

f{x) dx = TTife for k < n, and (A.l) 
— minimal for k > n. (A.2) 



j fix) Hk dx 

JR 



In (jA.2p the Hermite polynomials H^iy) = X!/=o ^ki are written in the vari- 
able y = (x — mi)/V2a rescaled by the width a = \/ rcii — m\. The truncated 
moment problem (jA.ip has to be augmented by the complementary condition 
(|A.2p for uniqueness. Functions solving the Eqs. (jA.ip are readily given by 



By virtue of the orthogonality of the Hermite polynomials, the first n coefficients 
Cfc = (iik) j — J fHkdx — X]f=o ^'s' directly determined by the known 

moments to/. If one was to truncate the series (jA.3|) after the n's coefficient, the 
resulting function /o may take negative values. If so, we use this negative part 
fo — fo 0(~/o) of the function /o = fo+fo to determine the higher coefficients 
to he Qk = — {Hk) J- for k > n. This yields a new function fi = fi + fi with 

a smaller negative part . Iteratively one approaches the desired solution f^o 
with arbitrary precision (cf. Fig. [6]). 



References 

[1] H.D. Zch, The Physical Basis of the Direction of Time 
(Springer, 2001). 

[2] D.N. Spergel, et.ai, ApJS, 148, 175 (2003); C.L. Bennett, 
Nature 440, 1126 (2006). 

[3] S.W. Hawking, Phys. Rev. D 32, 2489 (1985). 

[4] D.N. Page, Phys. Rev. D 32, 2496 (1985). 

[5] S.W. Hawking, R. Laflamme, and G.W. Lyons, Phys. Rev. 
D 47, 5342 (1993). 

[6] A.E. AUahverdyan and V.G. Gurzadyan, J. Phys. A 35, 
7243 (2002). 

[7] M. Castagnino, O. Lombardi, and L. Lara, Found. Phys., 
33, 877 (2003). 



25 



[8] G.N. Bochkov and Yu.E. Kuzovlev, Physica 106A 443 

(1981) and references therein. 

[9] D.J. Evans, E.G.D. Cohen, and G.P. Morriss, Phys. Rev. 
Lett. 71, 2401 (1993). 

[10] G. Gallavotti and E.G.D. Cohen, Phys. Rev. Lett. 74, 2694 
(1995); J. Stat. Phys. 80, 931 (1995). 

[11] J. Kurchan, J. Phys. A 31, 3719 (1998). 

[12] R. van Zon and E.G.D. Cohen, Phys. Rev. Lett. 91, 110601 
(2003). 

[13] U. Seifert, Phys. Rev. Lett. 95, 040602 (2005). 

[14] C. Jarzynski, Phys. Rev. Lett. 78, 2690 (1997); Phys. Rev. 
E 56, 5018 (1997); J. Stat. Phys. 98, 77 (2000). 

[15] L. Boltzmann, Ann. Phys. 57, 773 (1896). 

[16] W. Gibbs, Trans. Conn. Acad. 3, 229 (1875). 

[17] F. Debbasch, K. MalUck, and J.P. Rivet, J. Stat. Phys. 88, 
945 (1997). 

[18] J. Dunkel and P. Hanggi, Phys. Rev. E 71, 016124 (2005). 

[19] J. Dunkel and P. Hanggi, Phys. Rev. E 72, 036106 (2005). 

[20] A. Einstein, Ann. Phys. 17, 549 (1905). 

[21] A. Einstein, Ann. Phys. 17, 891 (1905). 

[22] D.J. Evans, D.J. Searles, Adv. Phys. 51, 1529 (2002). 

[23] C. Maes and K. Netocny, J. Stat. Phys. 110, 269 (2003). 

[24] G.E. Crooks, Phys. Rev. E 60, 2721 (1999); Phys. Rev. E 
61, 2361 (2000). D. Collin et al. Nature 437, 231 (2005). 

[25] H. Risken, The Fokker-Planck Equation (Springer, 1996). 

[26] C. Maes and M.H. van Wieren, Phys. Rev. Lett. 96, 240601 
(2006). 

[27] H. van Hees, V. Greco, and R. Rapp, Phys. Rev. C 73, 
034913 (2006). 

[28] Y.P. Raizer, Gas Discharge Physics (Springer, 1991). 

[29] F. Jiittner, Ann. Phys. 34, 856 (1911). 



26 



[30] S. KuUback and R. A. Lcibler, Ann. Math. Stat. 22, 79, 
(1951); S. KuUback, Information Theory And Statistics 
(Dover Publications, 1968). 



[31] J. A. Peacock, Cosmological Physics (Cambridge University 
Press, 1999). 

[32] N.D. BirrcU and P.C. Davics. Quantum Fields in Curved 
Space, 59 (Cambridge, 1994). 

[33] S. Weinberg, Gravitation and Cosmology, 508 (Wiley, 
1972). 

[34] M. Rigotti and F. Debbasch, J. Math. Phys. 46, 103303 
(2005). 



27 



